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Abstract 

It is showed that the very recently introduced Lagrangian loop 
formulation of the lattice Maxwell theory is equivalent to the Villain 
form in 2+1 dimensions. A transparent description of the classical 
loop action is given in pure geometrical terms for the 2 + 1 and 3 + 1 
dimensional cases. 



The loop formalism was introduced some time ago as a general analytical 
Hamiltonian approach based on the properties of the group of loops It 
provides a common framework for quantum gauge theories - it works for sev- 
eral models, as Maxwell theory [0], Chern-Simons theory 0, etc - and quan- 
tum gravity [Q. Furthermore, it works for Yang-Mills theories on a lattice P| 
and recently has been extended in such a way to take account of dynamical 
fermions 0. Very recently a loop action for the U(l) gauge theory has 
been built which lattice version leads to the Kogut-Susskind Hamiltonian. A 
Lagrangian formulation in terms of loops is interesting for multiple reasons. 
First, it offers the possibility of knit together the intrinsic advantages of the 
nonlocal loop description (the non-redundancy of gauge degrees of freedom 
and its geometrical transparency) and the computational power of numeri- 
cal simulations. This provides a a useful complement to the analytical loop 
studies. The classical action may be also relevant to perform semiclassical 
approximations. In this paper we show that the proposed action for the 
D=2+l dimensional case leads after a duality transformation to the discrete 
Gaussian model form. It is known that the same happens for the Villain 
form. Thus, for D=2+l, the minimal description provided by the loops ac- 
tually corresponds to that of Villain form instead of Wilson cosinus form. 
For D=3+l the loop action and Villain form give a similar description but 
the connection between them is more subtle. 

According to reference [0 the continuum Euclidean loop action for the 
abelian theory in D = d+ 1 space-time dimensions is given by 

S=^J dtd'x{-X^ai^)^X^ai^) + XMx)X^(x)} (1) 

where X^,{x.) is the loop current §fj (i?/*5^(x— y) and A is the three-dimensional 
Laplace operator. The time derivative of ^^(x) can be written as lim^^j^o 
It ^Ct+dtCt ^l/*^^(^~y)) where Ct denotes the loop Ct traversed in the opposite 
direction. 

In order to formulate (|1]) on a lattice we represent the continuum surface 
spanned by the loop C as a set of spatial loops Ct at different times t = 
0,1, ... ,T. We also replace the derivatives by finite difference operators and 
get 

t s ^ 
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where Xq^{s) is an integer which counts the number of times that the loop 
Ct traverses the hnk {s,i). 

We shall consider first the D=2 + 1 dimensional case. The loops Ct lie on 
a plane and have a unique decomposition in plaquettes Ct = YipP"''': where 
Up is the plaquette multiplicity, i.e. the number of times that plaquette p is 
contained in loop C. The first term in (|]) is then 



Y.-XUs)^Xlis) = Y.Y.-^P^P'X;{s)^X;,is) = Y.npnpd{p,p'^ (3) 

s p,p' s p^pl 

where D denotes the loop Ct+iCt, l{p,p') is the interaction between two 
plaquettes and is given by 

iip,p') ^Y.-xi>i-^)^x;,is) = E^;(^)G'(^,^')^M^') (4) 

s s,s' 

with 



/•27r ^2 iq{s-s') 

G{s, s') = / % _ r = 5,,,. (5) 

JO (27r)^ 2(2 — cosgi — cosg2) 

That gives 



Y: -X},is)^XUs) = Y: npnp,5p,y = Y.< = MD) (6) 

s p,p' P 

where A2{D) actually denotes the quadratic area, i.e. the sum of the squares 
of the constitutive plaquette multiplicities, of the surface S enclosed by the 
loop D. 

The second term on the action is 

E^c(5)^c(^) = A(C) (7) 

s 

where A(C) denotes the quadratic length -the sum of squares of link multi- 
plicities- of the loop C. 

We arrive then to the following geometric expression for the lattice action: 

s=^^Y.{MCt+iCt) + KCt)}. (8) 
^ t 
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Now, instead of expressing the loops as an integer field defined on the links 
of the lattice, we express them through the surface they enclose by giving an 
integer field defined on the spatial plaquettes. By duality this field is defined 
on the sites {x,y,t) dual to these spatial plaquettes so the configurations 
are given by integer fields k{x,y,t). The geometric terms in the action are 
expressed then as 



HCt) = E {(^(^ + 1' - KS:, y, t))' + {k{i, y + ^,t)- k{x, y, t)Y} (9) 

x,y,t 

and 

^2(^+1^) = E (^(*' y,t + ^)- y, t)Y- (lo) 

x,y,t 

Collecting these results, the action reads finally 

s = ^-^T.(Hs + f^)-Hs)? (11) 

which is the action of the integer Gaussian spin model. In this expression 
the covariance hidden in (|]) is manifest. 

On the other hand, the duality of the Villain form of the U(l) gauge 
theory to discrete Gaussian models -discrete Gaussian spin model for D=3, 
discrete Gaussian gauge model for D=4- was already established 0, |T0| 



Thus, for D=3, we have found that the minimal description that provides 
the loop action (2) is equivalent to the Villain form since that both lead by 
duality to (12). 

With the aim to explore this connection between loops and Villain form 
for an arbitrary dimension D we shall resort to the language of forms. Work- 
ing in this general framework we will show that we recover the previous 
equivalence for D=2+l and we will analyze in detail the relation between 
both theories at D=3+l. 

A p-form is a function defined on the p-cells of the lattice (p=0 sites, 
p=l links, p=2 cubes, etc.) over an abelian group which shall be R, Z, or 
U(l)=reals module 2tt. Integer forms can be considered geometrical objects 
on the lattice. For instance, a 1-form is a path and the integer value on 
a link is the number of times that the path traverses this link. V is the 
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co-border operator which maps p-forms onto (p+l)-forms. It is the gradient 
operator when acting on scalar functions (0-forms) and it is the rotational 
on vector functions (1-forms). We shall consider the scalar product of p- 
forms defined < a \ (3 >= J2cp C({c)(3{c) where the sum runs over the p-ccUs 
of the lattice. Under this product the V operator is adjoint to the border 
operator d which maps p-forms onto (p-l)-forms and which corresponds to 
minus times the usual divergence operator. The operator □ — Vd + dV is 
called the Laplacian and differs only by a minus sign of the current Laplacian 
A^A^. We shall use also the duality transformation which maps biyectively 
p-forms over (D-p)-forms. Under duality the border and co-border operators 
interchange. 

In the loop theory exposed before, physical configurations correspond 
to distributions of spatial loops at different times. These are spatial closed 
integer 1-forms. Then we consider only forms c which are for temporal links 
and closed [dc = 0). The border and co-border operators can be restricted 
to the spatial sections (t=constant) and we denote the spatial Laplacian as 
A. The loop action (2) can be expressed as 



-5 = y Z](< Ct+lQ, ^Q+lCt > + < Ct, Ct >) = y < c, — c> (12) 

and the path integral 

Z= Yl exp(-^<c,Hc>) (13) 
c 

(dc = 0) 

where c denotes a spatial loop i.e. c belongs to a t=constant slice. 

It is possible to solve the constraint 5c = by taking c — ds where s are 
integer spatial 2-forms. 

Z^j:'eM-^^<ds,^ds>) (14) 

The prime on the sum indicates that we must restrict the s's over which 
we sum since we have introduced the symmetry s — > s + dv ior arbitrary 
3-forms v (in other words, there are infinite surfaces with the same loop as 
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border). This restriction is the usual gauge- fixing problem and is solved more 
simply if we make a duality transformation in the spatial sections. In terms 
of the spatial (D-3)-forms *s (14) is written as 

Z^Yl e^P(-y < V * s, -^v * s >) (15) 

For D=2+l the dual of *s are 0- forms *s. We do not need to fix the 
gauge since a spatial loop in the border of a unique surface. 

If we use the Poisson sumation formula J2j f{j) — J2k I^oo d(j)f{(j))e^'^'^'^^ 
we can substitute the integer *s variables by integer 0-forms *m and real 

0- forms 0, and performing the Gaussian integration over we get 

27r2 1 

Z = ^exp( — < *m, — * m >) (16) 

m 9 

The above expression is a sum over point monopoles which interact through 
a Coulomb potential. Its form is equal to the one obtained from the Villain 
form of the U(l) theory. Those monopoles are defined in the centers of the 
spatial plaquettes and these configurations are in correspondence one-to-one 
with scalar forms on the three-dimensional lattice. In this case, as we have 
seen before, the loop theory describes exactly the topological excitation of 
the Villain theory. 

For the D=3-|-l case the dual of the s's in (14) are spatial 1-forms *s. 
To fix *s we impose *s{l) — for every link I & T, T being a maximal 
spatial tree. This selects only one *s among all *s's with the same value of 
V * s. If we use the Poisson sumation formula we can substitute the integer 
*s variables on the complementary of T, T by integer 1-forms *m and real 

1- forms (j). Once again the integration gives 

27r2 1 

Z = constant ^ exp( — < *m, — * m >) (17) 

*m ^ 
{d*m — 0) 

i.e. a sum over spatial closed loops which interact through a Coulomb 
potential. The Villain form of the U(l) gauge theory leads to the same 
expression but the monopole loops are not restricted to lie at t=constant. 
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Thus, the loop theory is not fully equivalent to the Villain U(l) theory at 
D=4. 

We can deepen our insight on the relation between the loop and the 
Villain theory following the opposite way: starting from Villain form and 
then repeating the previous treatment in terms of p-forms. The path integral 
for the Villain U(l) theory is 

Z = |(rf^)^exp(-||| V^-27rj|r) (18) 

j 

where 6* is a real periodic 1-form, that is, a real number 9i e [0, 27r] 
defined in each link of the lattice, j are integer 2-forms, defined at the lattice 
plaquettes, and || . |P=< •,■ >■ Using the Poisson sumation formula and 
integrating first in the continuum (f) variables and then over the compact 9, 
we finally get 



Z = (27r/?)-^^/^ E exp(-^<.,.>)cxi:exp(-^^(5)). (19) 
{ds = 0) 

where Np in the number of plaquettes of the lattice. We can now analyze 
the geometrical meaning of this equation, s are integer 2-forms, that is they 
are defined over bidimensional surfaces on the lattice. The condition ds = 
means to sum only over closed surfaces, the action < s, s > is proportional 
to the quadratic area of the surface. Selecting one direction of the lattice as 
a 'time' direction one can express < s,s > as (see Appendix) 



< s, s >= ^(< Ct,Ct> + < ct+ict, ^ct+ict > + < Vst, ^Vst >) (20) 

where represents a pure spatial 2-form. The first two terms in (20) 
correspond exactly to the loop action (12). If we work at D=2-|-l we have 
only one solution to dst = Ct+iCt- In this case the spatial co-border of a 2- 
form is 0, then the third term in (20) vanishes and we obtain the loop theory 
(12). 
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For D=3+l the loop action can not take account of the term < Vs^, -^Vs^ > 
The Lagrangian lattice loop description seems to miss some degrees of free- 
dom present in the gauge U(l) theory (the general monopole loops which are 
not contained in the t=constant hypersurfaces) . As it was previously men- 
tioned it has been shown that the loop action leads to the Kogut-Susskind 
Hamiltonian so for continuous time both descriptions are equivalent. 

It seems interesting to have a Lagrangian lattice loop description to con- 
tinue exploring some interesting points. For example, the pure 3-1-1 U(l) 
gauge theory was studied by means of Hamiltonian analytical method in the 
loop representations In this reference we found a continuous phase tran- 
sition for this model instead of the first order measured by means of standard 
Monte Carlo simulations. A possible explanation for this difference relies on 
the fact that we have used a cluster method which implies a formally infinite 
lattice which is free from the spurious topological effects occurring in finite 
lattices with periodic boundary conditions [0 |13[. This periodicity leads, 
on the confining side, to large monopole loops which wrap around the lat- 
tice. Those loops clearly have a fairly large action associated with them, so a 
jump is observed when they disappear at deconfining transition. This jump 
no longer survives in the infinite volume limit. The effect of this wrapping 
loops in driving the transition is crucial. Recently, Stack and Wensley [|14 



have calculated the heavy quark potential from the magnetic current due the 
monopoles. They have resolved the magnetic current into large loops which 
wrap around the lattice and simple loops which do not. They found that the 
long range part of the heavy quark potential can be calculated solely from 
the large, wrapping loops of magnetic current. 

The loop's method seems to be in a good position because the finite lattice 
is not one of its intrinsic features. The development of different Lagrangian 
computational techniques for studying the D=4 U(l) model in order to com- 
plement our previous Hamiltonian results is one of our present aims. Also, 
we are interested in extending the loop's action in such a way to include 
matter fields. 

We wish to thank, D.Armand-Ugon, R.Gambini, M.Grady and L.Setaro 
for useful discussions and comments. 

APPENDIX 

Selecting a direction on the lattice which we call time and consider the 
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division of the lattice in space sections of dimension D-1. Any surface 2-form 
s can be written s = J2t{^t + St) where s is purely spatial and s is temporal 
(its plaquettes have one link in the t direction joining the space sections t 
and t + 1). The t subindex correspond to each spatial slice of the lattice. If 
we consider the borders of these surfaces we have that they are always spatial 
1-forms (loops) in the slice t since the temporal parts of dsf cancel, dsf has 
two spatial components in the sections t and t + 1, {dst)inf and {dst)sup- We 
define q = {dstjsup- Then {dst)inf — —Ct- Imposing s to be closed we have 
dst = Ct+i — Ct- We note also that if p is a plaquette defined by a spatial 
link / and a temporal link, then the integers St{p) and Ct{l) coincide. The 
quadratic area is then 

< S, S >= ^(< St, St> + < St, St >) = 

t 

^(< St, St> + < dst, ^dst > + < Vst, ^Vst >) = 



^(< Ct,Ct> + < ct+ict, -^ct+ict > + < Vst, -Vst >) (21) 

In the second equality we have inserted in the second term DD^^ and 
used the adjoint relation between border and coborder operators. In this 
expression the operators belong to the spatial hypersurface, so they 'live' in 
dimension D-1. In the path integral the sum over s is replaced for a sum 
over the distributions of spatial loops Ct and spatial surfaces St such that 
dst = Ct+\Cf 
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